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Abstract 

The decay rate of Kl — > 7r°77 has been measured with the NA48 detector at 
the CERN SPS. A total of 2558 Kl — > tt°jj candidates have been observed with 
a residual background of 3.2%. The branching ratio is determined to be (1.36 ± 
0.03( S f a f) ± 0.03( SJ/S ^ ± 0.03( raorm )) x 10~ 6 and the vector coupling constant a v = 
—0.46 ± 0.03( sta t) ± 0.04/ syst y This result suggests that the CP-violation effects are 
dominating in the Kl — > 7r°e + e~ decay. An upper limit for the — > 7r°77 decay 
rate in the two photon mass region m 77 < m n o is also given. 
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1 Introduction 



The measurement of the K L — > 7r°77 decay is useful to constrain the CP 
conserving amplitude of the decay Kl — > 7r°e + e~ via two photon exchange. 
Previous measurements of the decay [1-4] have been compared with calcula- 
tions performed in the framework of the Chiral Perturbation Theory (xPT), 
the effective theory of the Standard Model at low energy in the hadronic sec- 
tor. These predictions are best described in terms of the two amplitudes A and 
B (referring to angular momentum states J = and J = 2 of the two photons 
respectively) in the Lorentz invariant expression for the double differential 
decay rate [5]: 
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with rriK and px the kaon mass and momentum, fc 3 and the two photon 
momenta and y max the kinematic bound for the y variable, given by: 
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The leading 0(p A ) xPT calculation predicts B=0, in qualitative agreement 
with the experimental observation of a z spectrum peaked at high values, but 
underestimates the Kl — * 7r°77 branching ratio by about a factor of three 
[6]. At 0(p 6 ) the rate and the m 77 spectrum can be reproduced by adding 
a contribution from the VMD mechanism [5,7,8], via the coupling constant 
a v [9-11]. The VMD mechanism could enhance the state J = 2 for the two 
photons, hence allowing a sizeable CP conserving contribution which is not 
helicity suppressed for the K L — > n°e + e~ decay [12]. 

In this paper we present a new measurement of the branching ratio and de- 
cay spectrum of the decay K L — > 7r°77, based on data collected in the years 
1998 and 1999 with the apparatus of the NA48 experiment. Quantitative in- 
formation on the VMD contribution is obtained from the spectrum at low z 
values. 



2 Experimental set-up 



The K L — > 7r°77 events are collected using the NA48 detector located at the 
CERN-SPS and primarily devoted to measure the parameter e'/e which char- 
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acterises direct CP violation in neutral Kaon decays into two pions [15] [16]. 
Neutral Kaons are produced in interactions of 1.4 x 10 12 450 GeV/c protons 
on a beryllium target during 2.4 seconds every 14.4 seconds. A system of 
sweeping magnets and collimators defines a neutral beam of 2 x 10 7 K L per 
burst and divergence of ±0.15 mrad. A fraction of the non interacting protons 
are redirected by channelling in a bent crystal, to a second target, 120 metres 
downstream of the first, to generate a Kg beam [17]. In order to tag Kg decays, 
the protons before hitting this target are detected by an array of scintillation 
counters (tagger). The two beams converge with an angle of 0.6 mrad at the 
calorimeter, « 120 metres downstream of the final collimators. The decay re- 
gion is contained in an evacuated cylindrical vessel 89 metres long separated 
from the NA48 detector volume by a thin Kevlar window 0.3% of a radiation 
length (X ). The neutral beam traverses the detector inside a 16 cm diameter 
vacuum pipe. A scintillation counter (AKS) placed on the K s beam defines 
the upstream edge of the fiducial region for decays from that beam. 
The most important detector element for this analysis is the quasi- homogeneous 
liquid Krypton electro-magnetic calorimeter (LKr) structured in 13212 read- 
out tower cells 27 X Q deep [18]. The ionization signal from each of the cells 
is integrated, amplified, shaped and digitised by 10-bit FADCs at 40 MHz 
sampling frequency. The energy resolution can be parameterised as: 

<j(E)/E ~ (0.09 ± 0.01)/ E © (0.032 ± 0.002)/^ © (0.0042 ± 0.0005), 

with E in GeV. The spatial and time resolutions are better than 1.3 mm and 
300 ps, respectively, for photons with energy greater than 20 GeV. 
Seven ring shaped counters (AKL), consisting of plastic scintillator and iron 
converters, surround the decay and spectrometer region in order to veto events 
with photons outside the calorimeter acceptance. Their efficiency has been 
estimated to be around 95%. The charged particles are reconstructed by a 
magnetic spectrometer consisting of four drift chambers (DCH) and a mag- 
netic dipole. The space resolution for each projection is 90 fj,m and the average 
efficiency is better than 99.5% per plane. For this analysis the spectrometer is 
used to veto events with reconstructed tracks. 



3 Data taking 

Events of the K L — > 7r°77 decay channel and of the channel K L — > 7r°7r°, 
used for normalisation, are collected by the same trigger, since their final 
states appear identical in all detector elements. The trigger decision is based 
on quantities which are derived from the sums of the energy deposited in 
the LKr calorimeter in groups of cells corresponding to 64 horizontal and 64 
vertical slices [19]. The trigger requirements are: 
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1) at maximum 5 peaks in each projection; 

2) total energy larger than 50 GeV; 

3) centre of gravity of the event, computed from the first moments of the energy 
peaks in the projections to be within 15 cm from the beam axis; 

4) proper time of the decay, computed from the second moments, less than 
5 tk s from the beginning of the decay region. 

The trigger efficiency has been checked to be better than 99.9% [16] using a 
minimum bias sample triggered by a scintillating fibre hodoscope placed in 
the LKr volume at a depth of 9.5 X . 

During the experimental runs roughly 180 Terabytes of raw data were recorded 
and pre-selected by an on-line software filter, whose criteria are subsequently 
tightened in the off-line analysis. 



4 Analysis 

The first step is to reconstruct the showers in the liquid krypton calorimeter 
summing the energies of the cells in a circle of 11 cm radius from a local 
maximum. The time and the position of the shower are derived from the 
most energetic cells and from the centre of gravity of the energy. A partial 
overlap between two showers is resolved using the shower profile obtained 
from simulation. The shower energy is corrected for: 

i) the residual 0.4% non- uniformity in the calibration electronics chain, com- 
paring energy and momentum of electrons from K e3 events as a function of 
the hit cell. The uniformity is improved to the 0.15% level; 

ii) the 1% peak-to-peak modulation of the response as a function of the impact 
point in the LKr cell due to the effects of a 2 mm gap between electrodes 
in the vertical direction and to the finite integration time in the horizontal 
one; 

iii) the energy-leakage in the beam pipe at small radius, the material before the 
active volume and the presence of not properly working cells; 

iv) the small (0.1%) non-linearity. 

Since the Kl — > 7r°77 event is characterised by four clusters with the only 
constraint of two showers coming from a it decay, all the events with at 
least four clusters are considered. Each cluster must be located well inside the 
calorimeter volume, i.e. within an octagon of 113 cm apothem, with 15 cm 
inner radius. It must be at a distance larger than 2 cm from any non working 
cell and separated from another cluster by at least 10 cm. The cluster energy 
must be between 3 GeV and 100 GeV and the time within 3 ns of the event 
time, defined as the average of the cluster times. The total energy of the 
clusters must be in the 70-170 GeV interval. 
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The two photon invariant mass is: 



m wi = ~ — ( 2 ) 

ZLKr — Zk 



where Ei, Ej are the cluster energies, dfj = [(xj — Xj) 2 + (y, — Vj) 2 ] with Xi 
and Hi the horizontal and vertical positions of the i-th cluster at the LKr front 
face. The longitudinal coordinate of the decay vertex, zk, is computed from 
the electro-magnetic calorimeter information assuming the nominal kaon mass 
for the parent particle: 



Z K — ZLKr , [o) 

m K 



where zlkt is the longitudinal coordinate of the LKr front face. The overall 
energy scale is fixed by the fit of the position of AKS counter, that vetoes all 
Kg decays occurring upstream [15]. 

We only consider events for which there is no proton detected in the K$ tag- 
ger in time (within ±2 ns) with the event; this rejects 10.65% [16] of the Kl 
decays, due to accidental coincidences with a proton in the tagger. 
Because of the small branching ratio, the K L — > 7r°77 signal is subject to 
large backgrounds such as Kl — > 7r 7r°7r° decays, with missing and/or over- 
lapping photons, mis-reconstructed Kl — > tt tt° decays with 7 conversion in 
the spectrometer and events generated by accidental pile up of particles from 
two different events. We distinguish two classes of cuts in the analysis: 2n° 
rejection and 3n° rejection. 



4-1 2ir° rejection. 

Three different combinations of 77 pairs can be formed from four showers. A 
X 2 variable is defined: 



x 2 



(mi >2 + m 3A )/2 - m^o 



+ 



(mi j2 - m 3A )/2 



(4) 



where a± are the resolutions of (m^ ± m 3A )/2 measured from the data 
and parametrised as a function of the lowest photon energy, typically a + ~ 
0.42 MeV/c 2 and a_ w 0.83 MeV/c 2 . In order to select a 7r°7r° candidate the 
best combination is required to have x 2 ^ 13.5, while the 7r°77 events are 
selected with x 2 > 300. We retain for further analysis all combinations that 
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give a photon pair (that we call m^) within 3 MeV/c 2 from m n o and the 
other pair (corresponding to 777.3,4) outside the window 110 — 160 MeV/c 2 . 
The 27T° background concentrates in the most interesting region for the de- 
termination of a v , namely at low m 3j4 value (m 3i4 < 240 MeV/c 2 ) where the 
Kl — > 7T°77 signal is smaller. Due to mis-reconstructions one of the 7r° may be 
reconstructed with a wrong mass. This can happen in case of a 7-conversion or 
7T° Dalitz decay with one electron outside the calorimeter acceptance that bi- 
ases the photon energy reconstruction. We use the first three drift chambers as 
a veto against conversions and Dalitz decays if there are hits in coincidence, 
±20 ns, with the event time. In addition we reject events with an overflow 
condition in the readout of more than one plane of the drift chambers. This 
condition is defined by the occurrence of more than seven hits in a DCH plane 
in any 100 ns time interval. In this case the front-end readout buffers are reset 
and the DCH information cannot be used reliably. 

The remaining 2tt° background is reduced requiring p T 4 > 40 MeV/c, where 
Pt,4 is the transverse momentum to the kaon direction of 74, which is the 
lowest-energy photon of the non-7r° pair of photons. The variable pt,a is on 
average smaller for 27r° events than for 7r°77 events and even lower when the 
photon energy is underestimated. Although this requirement reduces the ac- 
ceptance for 7r°77 at high y and low m 3 4 , it does not degrade the sensitivity 
to the B amplitude as the (y 2 — ?/^ ax ) factor in equation 1 is already small. 



4-2 37T° rejection. 



In order to reduce the background from Kl — > 7r 7r°7r° decays, we require 
the event to have only four clusters in-time, within 3 ns, and no activity in 
the AKL veto counters in coincidence with the event. As most of the 3ir° 
background events have one or two photons outside the LKr acceptance, the 
reconstructed zk according to formula (3) moves downstream the beam axis 
by at least 6 metres. Therefore, the fiducial volume is restricted to the first 
30 metres of the K L decay region. Events with missing energy are also reduced 
by requiring the center of gravity, computed from the first momenta of the 
clusters (Rcog), to be less than 4 cm from the beam axis. 
The simulation of ~ 3 x 10 9 Kl — > 7r 7r°7r° events shows that after applying 
these cuts the background is composed of events with missed photons: 36% 
of the sample has two missing photons, while 61% of the sample results from 
one photon missing the detector and two showers overlapping in the LKr. The 
remaining events have all the six photons inside the calorimeter but only four 
reconstructed because of two overlaps. The corresponding showers are broader 
than for a single isolated photon. 

Events with missed photons are rejected on the basis of a variable, Z max , 
designed to estimate, for these events, the true kaon decay vertex position. 
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The same estimate, when applied to signal events, produces an unphysical 
result, with a distribution showing many events in the region upstream of 
the final collimator. The distribution of this variable is shown in figure 1. 
A condition Z max < — 5 m rejects almost 99% of the 3tt° background while 
retaining 46% of the signal. 

The procedure to define Z max is based on the analysis of the possible topologies 
of events with missed photons, as shown in the following figure: 
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cases 1, 2 and 3 : 

a 7T° can be reconstructed and the true vertex coordinate z 77 is obtained 
using equation (2) with m 7i7j = m n o. There are five possible 77 combi- 
nations, other than the one of the reconstructed ir°, but we consider only 
those that have z 77 smaller than {z K — 6m). We call Z % the vertex estimate 
which is closest to this limit. 
case 4 : 

one 7T° has a missing photon and two photons from the other 7r°'s overlap. 
Starting from the hypothesys of two 7r°s which share a photon, we compute, 
for each group of three photons, the vertex coordinate £ 777 , assuming each 
photon in turn to be the superposition of two photons with the same position 
in the calorimeter. For example if photons 1, 2 and 3 come from two 7r°s 
and photon 2 is overlapped, then ^ 777 is computed as: 



1/1 1 | 

hl1 ~ ZLKr ~ ^ \E^Ai + E^dlJ (5) 

Out of the twelve possible combinations we consider only the triplets with 
an invariant mass below (m K — m^o) and the vertex estimate closest to, but 
smaller than, [z K — 6m). 

For each event we define the variable Z max = max(Z 11 , Z 777 ). The residual 
3n° background originates from two overlaps of photons such that the Zk is 
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Fig. 1. Distribution of the variable Z max used to reject the 3ir° background with 
missing photons, as it looks after all the other requirements are applied. Dots are 
7r°77 candidates and solid line Monte Carlo simulation. The arrow indicates the 
region. 

reconstructed unbiased, but the shape of one or two showers is different from 
that of a single isolated photon. As estimator of the cluster width we use the 
RMS of its energy profile in x and y: 



W x{ y) 



i G 5 x 5 cells around the impact cell (6) 



where Ax(y)t is the distance in the x(y) direction of the i-th cell from the 
center of gravity of the cluster energy. An energy- dependent cut is applied 
to the width of each cluster. The value of the cut is defined in terms of the 
average w and the standard deviation a of the distribution of w measured in 
photons from good 27r°. We require w < (w + 1.8a) on the photons of m 3i4 
and w < (w + 3a) on the photons of m^. The signal reduction is at the level 
of 15%. 



5 Background evaluation 



The Kl — > 7r°77 candidates that survive all the above selections are 2558 in 
the signal region defined as 132 < m 12 < 138 MeV/c 2 . For 345 of these events 
more than one combination satisfies the 7711,2 cut and therefore the mass as- 
signment is ambiguous. 

We estimate the contamination from 27r° decays and pile-up processes directly 
from the data. To quantify the — ► 7r°7r 7r° background we rely on Monte 
Carlo. 

We study the Kl — > 7r°7r° background in the signal region using K$ — > 7r°7r° 
decays from different Ks data samples: the proton tagged Ks sample, col- 
lected concurrently with the Kl data, and a special Ks high intensity run 
(K S HI) taken in 1999. In the case of the K s tagged sample we ask for a 
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tight coincidence within ±0.5 ns between the tagger and the LKr signals. 
In this way we enrich the number of 2tt° events contributing to the back- 
ground. Conversely the only contribution to 7r°77 in this sample comes from 
the Kl contamination due to accidental coincidences in the tagger, which is 
estimated to be 3.5% with the cut used. After applying the 7r°77 selection 
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Fig. 2. 7713,4 distribution for the 3.5% of the Kl data (dots) and for the candidates 
found in the Ks tagged sample (histogram). No excess is visible in the low-mass 
region 

criteria, the 7773,4 invariant mass distribution is shown in figure 2 for events 
in the tagged sample compared to the original Kl — ► 7r°77 mass distribution 
scaled by the Kl contamination factor. The 2n° background is estimated from 
the tiny excess of tagged events measured where the 7r°77 signal is smaller, 
namely for 7773,4 < 0.240 MeV/c 2 while for 7773,4 > 0.240 MeV/c 2 it is lin- 
early estrapolated from control regions on the sidebands of the ir° peak ( 
127 < 777,1,2 < 130 MeV/c 2 or 140 < m 1>2 < 143 MeV/c 2 ). This procedure 
leads to a background evaluation of (4.1 ± 1.6( s tat) ± lA^gyst)) events. Here the 
systematic error is related to the subtraction of the Kl contamination and to 
the fact that K s events have a distribution of the decay vertex populated only 
in the region between 6 and 18 metres. 

The high intensity Ks data sample provides about 60 % of the tagged Ks — > 
7T 7r° statistic with negligible Kl component. The number of background 
events found is in agreement with the previous estimate. 
The cuts applied to reduce the Kl — > 7r 7r°7T background are mainly of ge- 
ometrical nature and can safely be simulated. The most critical cut is the 
one on the cluster width, for which an accurate description of the LKr re- 
sponse for overlapping showers is needed. For this last variable we compare 
data and Monte Carlo distributions in figure 3. The simulation describes the 
data within 3% both in shape and quantity being normalised to the Kaon flux. 
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5 10 15 20 2: 
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Fig. 3. Distribution of the w variable used to reject the Sir background (shaded 
area) with overlapped showers after all the other requirements. The white histogram 
corresponds to the sum of signal and background distributions. Monte Carlo data 
(solid line histogram) follow the experimental data (dots). 

After the cut the estimated residual background is (70.2 ± 8.2( stat ) ± 7.0( syst )) 
events, where the systematic uncertainty comes from the non-Gaussian tail of 
the calorimeter response mainly due to photo-hadron reactions. 
Finally, the pile-up events, which could originate for example from two over- 
lapped in-time kaon decays, can be quantified comparing the normalised dis- 
tributions of the center of gravity at the LKr for Kl — > 7r°77 candidates 
and for a good 27r° data sample. An excess at Rcog > 4 cm is found in the 
Kl — > 7r°77 distribution, after having subtracted the 37r° background (figure 
4). When the excess is linearly extrapolated under the Kl — > 7r°77 signal 
we obtain a pile-up background level of (8.1 ± 3.5( sta t) ± 4.0( SJ/St )) events. The 
choice of extrapolation function is dominating the systematic error; however 
the linear one is expected from a uniform center of gravity distribution and is 
confirmed by the Rcog distribution for events with out-of-time clusters. 

In Table 1 is reported the summary of the background evaluation in the signal 
region and, as cross-check, in the control regions defined above. 

The overall background is 3.2%. The invariant mass distributions mi t 2 and 



sample 


7T°77 


3tT° 


2vr° 


Pile-up 


Total bck % 


signal 


2558 


70.2±10.7 


4.1±2.1 


8.1±5.3 


82.4±12.1 3.2 


sidebands 


44 


30.6± 5.9 


2.1±1.6 


4.8±3.5 


37.5± 7.0 



Table 1 

Number of candidate and background events in the signal region (132 < m\i < 
138 MeV/c 2 ) and in control regions (127 < mi 2 < 130 MeV/c 2 or 140 < mi 2 < 
143 MeV/c 2 ) on the sidebands of the ir° peak. 

m 3j4 are shown in figure 5. The background estimate is in agreement with that 
found in the two control regions where 44 background events are observed in 
the data with an expectation of 38 ± 7. 
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Fig. 4. Rcog distribution (left) for the candidate events (dots) for 2tt° events (open 
circle). The residual background in the tail of high RcOG is linearly extrapolated 
below the peak as suggested from the Rcog distribution of events with out-of-time 
clusters (right). 
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Fig. 5. mi 5 2 and mass distributions for the Ki — > ir°jj candidates events (dots) 
and Monte Carlo data (solid histogram) where the fitted value of a v has been used in 
the event generation. The shaded histogram is the sum of the three background types. 
In case of ambiguous mass assignment, only the combination with mi. 2 closest to 
m_o mass is included in these distributions. 



6 Results and discussion 



As confirmed by the Monte Carlo, most of the events in the low mass 
region shown in figure 5 are genuine Ki — > 7r°77 events which suffer from the 
wrong combination of photons. Events with an ambiguous assignment for m 12 
have a flat distribution in the 7713,4 mass region which will inflate the amount 
of events in the low-mass m 3j 4 < 0.240 Mev/c 2 region. To extract a value 
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Table 2 

7r°77 non- ambiguous candidates, number of background events and acceptance in 
20 bins of m^^ in the range of [0 — 400] MeV/c 2 . The normalisation factor, i.e. 
number of decays in the fiducial volume N(Kl), is equal to 23.9 x 10 9 . This number 
is based on the number of measured N(Kl — ► 7r°7r° ) decays. 



for the vector coupling constant a v we exclude events with ambiguous mass 
assignment. The events in the low- mass region are thus reduced by ~ 50%. In 
Table 2 we report the numbers of 7r°77 candidate, background and acceptance 
in the 20 bins of the [0 — 400] MeV/c 2 m 3 4 interval. We perform a fit of the 
bidimensional distribution of the two relevant kinematic variables m 3i4 and 
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y — \E 3 — E^l/niK, defining a likelihood function as follows: 

InL = [ N i(y, m n) ' inEi(y, m 77 ) - E^y, m 77 )] (7) 

where Ni(y,m iy ) is the measured number of events in the i-th bin while 
Eiiy^m^) is the corresponding expectation value given by the simulated 
Kl — > 7r°77 plus the estimated background. The model adopted for the sim- 
ulation of the Kl — > 7r°77 process is described in [9] and it is normalised to 
the total number of Kl decays in the fiducial volume. The uncertainty of the 
parameterisation of the Kl — > 3n vertex [13] is included in the systematic 
error. The best fit, with a x 2 of 31.1 for 30 degrees of freedom, gives: 

a v = -0.46 ± 0m {stat) ± 0.0A {syst) 

The systematic error includes the following contributions: 0.02 from the un- 
certainty in the parameters of the theoretical model, 0.03 from the uncertainty 
of the acceptance evaluation and 0.02 from background evaluation and sub- 
traction. 

In figure 6 we show the mi^ distribution together with the evaluated back- 
ground and the 0(p 6 ) xPT expectation for the fitted value of the vector 
meson exchange parameter a v in the three low-mass regions a) m 3t 4 G [30 — 
110] MeV/c 2 , b) m 3i4 G [160-240] MeV/c 2 and c) m 3A G [240-260] MeV/c 2 . 
Our data clearly disfavour a v = 0. 

The branching ratio is obtained from the number of 7r°77 candidates, includ- 
ing those with ambiguous mass assignment, normalising to the number of 7r°7r° 
events observed in the same sample of K L decays and selected by the same 
cuts, in order to minimise the uncertainty in the acceptance corrections and 
is found to be: 

BR(K L -> 7T°77) = (1.36 ± 0.03 {stat) ± 0.03 {syst) ± 0.03 (norm) ) x lO" 6 

Both the correction due to the Dalitz decay difference between the 27r° normal- 
isation mode and the 7r°77 decay and the one due to the K L —K S interference in 
7r°7r° decays [16] are taken into account. The uncertainty related to the experi- 
mental knowledge of the K L — > 7r°7r° branching ratio [21] is quoted separately. 
The computed average value of the acceptance for the 7r°77 process is 7.62% 
and for the 7r°7r° is 8.72%. The uncertainty related to this computation is the 
main source of the 2.5% systematic error in the branching ratio measurement 
since, as shown in Table 2, it has sharp variations. The residual background 
estimation, the uncertainty in a v and the exclusion of the ambiguous events 
each gives less than 1% variation on the branching ratio value. 
As shown in figure 6 there is no evidence for a signal in the invariant mass 
region m 3i4 < m n o where we count 9 7r°77 candidates and (4.5 ± 2.7) back- 
ground events. Using events with y < 0.2 in 5 bins of m 3 4 between 30 and 
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Fig. 6. Comparison between data (crosses) and the 0(p G ) expectations for 
a v = —0.46 and a v = (solid line histogram) in the mi 2 variable for the three 
regions of the low-mass spectrum: a) 7713,4 G [30 — 110], b) m^^ £ [160 — 240] and 
c) w<3.4 G [240 — 260] . Background is shown in the shaded area histogram. 

110 MeV/c 2 , we compute a model-independent upper limit [20]: 

BR(K L -> vr 77)| m3i4e[ 3o__iio] Afe»/©», w e[o-o.2] < 0.6 x 10~ 8 90% C.L. 

For this calculation we use the acceptance computed from events generated 
according to phase space, in a restricted region where it is almost uniform. 
Given the negligible contribution of the amplitude A(J = 0), our measurement 
of the vector coupling constant a v allows the CP-conserving component of the 
K L — > 7r°e + e _ decay [12] [14] to be predicted as: 

BR{K L -> 7r e + e-)\cp-conservin 9 = (4.7±f;l) x 10- 13 

This suggests that the CP-violation effects dominate in the Kl — * 7i°e + e~ 
mechanism. 
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